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Wc apply the 3-dimcnsional radiative transport eode Wind3D to 3D hydrodynamic models of 
Corotating Interaction Regions to fit the detailed variability of Discrete Absorption Components 
observed in Si IV UV resonance lines of HD 64760 (BO. 5 lb). We discuss important effects of the 
hydrodynamic input parameters on these large-scale equatorial wind structures that determine 
the detailed morphology of the DACs computed with 3D transfer. The best fit model reveals 
that the CIR in HD 64760 is produced by a source at the base of the wind that lags behind 
the stellar surface rotation. The non-corotating coherent wind structure is an extended density 
wave produced by a local increase of only 0.6% in the smooth symmetric wind mass-loss rate. 





Figure 1: Schematic representation of two CIRs in 
the plane of the equator used in 3D radia- 
tive transfer calculations with Wind3D. 



stellar radii. The CIRs can be produced by inten- 
sity irregularities at the stellar surface, such as dark 
and bright spots, magnetic loops and fields, or non- 
radial pulsations. The surface intensity variations 
alter the radiative wind acceleration locally, which 
creates streams of faster and slower wind material. 

We investigate to what extent the CIR wind model 
can quantitatively explain the detailed DAC prop- 
erties observed in BO. 5 lb supergiant HD 64760. 
Fullerton et al. 1997 pointed out how exceptional 
lUE data from the MEGA campaign (Massa et 
al. 1995) have made it a key object for study- 
ing the origin and nature of hot-star wind variabil- 
ity. We discuss the development of a new 3D ra- 
diative transfer code Wind3D to best fit the DAC 
morphology (e.g. shape and flux changes) observed 
in the unsaturated absorption portion of the Si IV 
A1395 line. The hydrodynamic CIR models devel- 
oped for this purpose are discussed in a compan- 
ion paper (Blomme 2007) . Animations are online at 
alobel . f reeshell . org/conf erence . html. 



1 Introduction 

Discrete Absorption Components (DACs) observed 
in the broad P Cygni profiles of UV resonance lines 
are important tracers of the dynamics of line driven 
winds in massive hot stars. DACs are observed to 
propagate bluewards through the UV line profiles on 
time scales comparable to the stellar rotation period 
(Massa et al. 1995; Prinja 1998). Hydrodynamic 
models of Corotating Interaction Regions (CIRs) 
have been proposed by Cranmer & Owocki (1996) to 
explain the observed DAC properties qualitatively. 
These large-scale wind structures are spiral-shaped 
density- and velocity-perturbations winding up in 
or above the plane of the equator that can extend 
from the stellar surface to possibly several tens of 



2 3D radiative transfer: WindSD 

Wind3D computes 3D spatial non-LTE radiative 
transfer in the 2-level atom approximation for opti- 
cally thick resonance lines formed in scattering dom- 
inated extended winds of massive hot stars. Its im- 
plementation is based on the finite element method 
described by Adam (1990). The code accepts arbi- 
trary 3D wind-density and -velocity distributions. 
The 3D transfer scheme further extends Adam's 
Cartesian method with three new aspects: (i) We 
considerably accelerate the (exact) lambda iteration 
of the source function on 71^^ grid points with appro- 
priate starting values from the Sobolev approxima- 
tion, (ii) Since the lambda iteration is the bottleneck 
of the numerical transfer problem we fully parallelize 
the mean intensity integration over 80^ spatial an- 
gles, (iii) We introduce a new technique that 3D in- 
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Figure 2: Hydrodynamic CIR models for spot opening angles $sp=:90° (left-hand panel) & 180° (right- 
hand panel). The dynamic spectra (middle panels) reveal width changes of the DAC base for an 
observer viewing the rotating models from the south side of the images in the equatorial plane. 



terpolates the converged (non-Sobolev) source func- 
tion to a higher resolution (spatial) grid of 700^ grid 
points to solve the final 3D transfer equation for very 
narrow line profile functions. This method resolves 
the small flux variations in the absorption portions 
of very broad unsaturated P Cygni profiles. 

Parameterized 3D models of CIRs already provide 
comprehensive comparisons to general properties ob- 
served in DACs. We consider a /3-velocity law for an 
isothermal wind with /3 ~ 1. The smooth wind is 
perturbed with 3D spiraling density enhancements 
wound around the central star (Fig. 1). The wind 
velocities inside the CIRs also assume the /?-law and 
are directed radially (velocity vectors drawn in the 
equatorial plane). The CIR model of Fig. 1 causes 
the width of DACs computed for an observer in the 
plane of the equator to decrease because the range 
of velocities in the CIR projected in the observer's 
line of sight (inside the cylinder in front of the stellar 
disk) decreases at larger distances from the star. 

For the hydrodynamic CIR models a local radia- 
tion force enhancement (or surface 'spot') is intro- 
duced at the base of the stellar wind (Blomme 2007). 
The resulting equatorial density- and velocity- 
structures are inserted around the star with a thick- 
ness of 1 i?* around the equatorial plane. Outside 
this region, the model density and velocity assume 
the smooth wind values. In Fig. 2 the spot opening 
angle <i>sp is increased from 90° to 180°. The spot 
CO- rotates with the stellar surface Wsp=Wrot, and the 
spot intensity Asp=0.5. The increase of $sp con- 
siderably alters the FWHM evolution of the DAC 
computed in Si iv A1395. The DAC base broadens 
because extra wind material injected by the spot be- 
comes more spread out over the equatorial plane. 



The maximum oi p/ pa in the CIR occurs within ^^5 
i?* above the stellar surface. Inside this region extra 
wind material is distributed over a larger geometric 
region above the bright spot, which also consider- 
ably broadens the density contrast in the tail of the 
CIR. When this region rotates in front of the stel- 
lar disk (around rotation phase 0.2 in the dynamic 
spectra) the DAC base broadens because the range 
of velocities projected in the observer's line of sight 
(that contribute to the DAC opacity) increases. The 
wind streams almost radially through the CIR den- 
sity structure. The wind flow decelerates due the rel- 
ative increase of wind density in a region behind the 
CIR density contrast maximum. It causes a trailing 
velocity plateau (between the star and CIR) which 
strongly contributes to the DAC absorption at larger 
distances from the star since the line optical depth 
in the Sobolev approximation is r ~ p / \dv/dr\. 

3 Modeling DACs in HD 64760 

We determine a recurrence period of 10.3±0.5 d for 
the two DACs (right-hand panel of Fig. 3) observed 
in the Si iv fines of HD 64760 in 1995 January 13- 
29 (Massa et al. 1995) by flux flltering the nearly 
horizontal rotational modulations (Lobel & Blomme 
2007). We can assume that this star is observed 
equator-on (sin i~l). The surface rotation velocity 
of the fast-rotating supergiant is Wrot=265 kms~^ 
which, for R^,=22 Rq, yields a rotation period of 
4.13 d. The latter period is 2.5 times shorter than 
the observed DAC recurrence period. It reveals that 
the spot cannot co-rotate with the stellar surface. 
We therefore use a single bright spot with the spot 
parameter Vgp set equal to Uiot / 2.5 in the hydro- 
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Figure 3: Best fit dynamic spectrum (left-hand panel) of Si IV A1395 compared to the observed spectrum 
(right-hand panel) of HD 64760. The computed DAC shape fits the observations in detail. 



dynamic wind models. We compute a large grid of 
dynamic spectra from a grid of hydrodynamic mod- 
els for a broad range of spot parameters 4>sp and 
Asp, using the smooth wind properties of HD 64760. 
We obtain a best fit from a least-squares fit proce- 
dure between the observed and computed DACs for 
$sp==50°±5° and ^p==0.1±0.05 (see the hydrody- 
namic model in Fig. 1 of Blomme 2007). In the 
left-hand panel of Fig. 3 the velocity positions of 
the flux minima in the computed DAC differ by less 
than ~50 kms~^ from the observed velocity posi- 
tions (white dots) for d < T < 3.5 d, and 10 d 
< T 15.5 d. The FWHM of the computed DAC de- 
creases from ~100 kms~^ at T=0 d to ~20 kms~^ 
around r=3.5 d, in agreement with the narrowing of 
the observed DAC. The DAC width remains almost 
constant over the following 6.5 d, after which it fades 
away. The 'tube-like' extension of the DAC base is 
also observed (right-hand panel). This characteris- 
tic DAC morphology can only correctly be computed 
with hydrodynamic structured wind models. 



4 Conclusions 

We demonstrate with 3D radiative transfer calcu- 
lations in hydrodynamic CIR models of HD 64760 
that the DACs observed in the Si iv UV resonance 
lines are due to a region of enhanced mass-loss at 
the base of the wind that lags 2.5 times behind the 
surface rotation. Our detailed DAC modeling re- 
veals that this region is not locked to the stellar sur- 
face. It indicates that magnetic fields at the equator 



are an unlikely source for additional wind material 
yielding an asymmetric structured wind with an ex- 
tended density wave in the equatorial plane. The 
integration of the best fit hydrodynamic CIR model 
for HD 64760 provides a very small increase of only 
0.6% in the mass-loss rate of the spherically sym- 
metric smooth wind model. It signals that DACs 
are generally expected in hot-star winds since they 
result from rather small variations of the spheri- 
cally symmetric mass-loss rate. These coherent CIR 
structures may become perturbed by wind clumping 
on much smaller length scales. They will however 
built up again very rapidly as well, provided that 
the perturbation time-scales are sufficiently short for 
the large-scale wind structures to completely develop 
and to produce the slowly blueward drifting recur- 
ring DACs in unsaturated UV spectral lines. 
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